Oscillator state reconstruction via tunable qubit coupling in Markovian environments 
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We show that a parametrically coupled qubit can be used to fully reconstruct the quantum state 
of a harmonic oscillator, even when both systems are subject to decoherence. By controlling the 
coupling strength of the qubit over time, the characteristic function of the oscillator at any phase 
space point can be directly measured by combining the expectation values of two Pauli operators. 
The effect of decoherence can be filtered out from the measured data, provided a sufficient number 
of experimental runs is performed. In situations where full state reconstruction is not practical 
or not necessary, the method can still be used to estimate low order moments of the mechanical 
quadratures. We also show that in the same framework it is possible to prepare superposition states 
of the oscillator. The model is very general but particularly appropriate for nanomechanical systems. 



INTRODUCTION 

With an increasing number of experiments bringing 
mesoscopic oscillators to the quantum regime, it has be- 
come essential to develop efficient readout methods that 
are able to verify quantum effects in such systems. The 
ultimate goal is the full reconstruction of the system's 
density matrix. Most of the current schemes to measure 
nanomechanical motion consist of weak continuous mea- 
surements of the oscillator's position (see e.g. [1]). Even 
though these techniques are useful for the monitoring of 
ground state fluctuations, they are not easily applica- 
ble to arbitrary nonstationary states, where the relation 
between the quantum state and the measured data is ex- 
tremely complicated. Recently a method that allows full 
reconstruction of an arbitrary state of a mechanical os- 
cillator has been proposed in the context of optomechan- 
ics [H, however it relies on tomography techniques that 
involve transformations of the measured data. In both 
cases it is challenging to identify an unknown quantum 
state with high accuracy. 

The idea that a two-level system (or qubit) could be 
used as a probe to reconstruct the quantum state of a 
continuous-variable system has been initially proposed 
in the field of cavity- QED, due to the impossibility to 
directly measure the intracavity field [sl, Q . Similar ideas 
have subsequently been introduced in the field of circuit- 
QED [H , ion traps [6] and nanomechanical oscillators . 
In all these works, information about the initial state of 
the oscillator is inferred from the transition probabilities 
of the qubit, after the two systems have interacted. Typ- 
ically, the value of the characteristic function |8| of the 
oscillator at several phase space points can be directly 
extracted from the qubit statistics. 

In this paper we report a further advance in this di- 
rection. Firstly, we assume to control the parametric 
coupling between the qubit and oscillator as a function 
of time, a very interesting possibility opened by some 
recent experimental and theoretical works (ol-llij. With 
this choice, any point of the oscillator's phase space can 
be reached by appropriately tailoring the time depen- 



dence of the coupling. The corresponding value of the 
characteristic function is then measured by combining 
the expectation values of two Pauli operators of the qubit. 
Differently from methods involving resonant interactions 
[sl, [HI, 0, no auxiliary degrees of freedom or additional 
displacement operations are required in our scheme. 

The second contribution we give to the field is the full 
inclusion of decoherence in the treatment, since to the 
best of our knowledge the effect of noise on such state 
reconstruction procedures has yet to be explored system- 
atically. The inclusion of decoherence is essential espe- 
cially when dealing with mesoscopic systems, where the 
coupling of both the qubit and the oscillator to a ther- 
mal environment cannot be neglected. When the envi- 
ronment is described via the standard Markovian master 
equation, we obtain the remarkable result that the char- 
acteristic function can still be reconstructed in full detail. 
The price to pay is that a larger number of experimen- 
tal runs will be needed to extract the same amount of 
information as compared to the ideal decoherence-free 
case. In practice, decoherence will limit the size of the 
phase space region where the characteristic function can 
be measured in a reasonable number of runs. When only 
the nearbies of the origin can be accessed, our method 
still gives useful information by providing low order mo- 
ments of the mechanical quadratures. 

As an addition we show, and examine, how well our 
system can be used to prepare superposition states of the 
oscillator in the presence of decoherence, and we provide 
a detailed analysis of the resulting characteristic function. 

Having proposed a method to reconstruct the charac- 
teristic function, we remark that it constitutes a very 
useful description of a quantum state, even though it is 
sometimes overlooked in favor of the Wigner represen- 
tation [8]. As we emphasize in this paper, a number of 
properties such as purity, squeezing or the presence of 
coherent superpositions can be investigated from it as 
directly as from the Wigner function. Since any moment 
of the mechanical quadratures can be extracted from it 
in a straightforward way, the characteristic function is 
the ideal representation to invesigate statistical proper- 
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ties such as squeezing or non-Gaussianity [1^. Moreover, 
an efficient method to estimate nonclassicality directly 
from a finite collection of measured characteristic func- 
tion values has been recently proposed [Tij . 



THE MODEL 

We consider a qubit with ground state |^), excited 
state |e), bare energy Uq and tunneling energy (5, coupled 
to a harmonic oscillator of frequency Q. The coupling 
strength g{t) is for now a generic function of time. The 
Hamiltonian, in units of = 1, is 



where po is the state of the harmonic oscillator that we 
want to reconstruct, and |+) = -^{\9) + [2o[. After 
an interaction time t the system evolves to a state ptot {t) 
according to (|3|). At this point, we measure either the 
observable ax or a^y, which completes a single "run" of 
the experiment. A sufficient number of runs will have 
to be performed in order to estimate the average values 
i^ji'^)) — "tr {ptot(^)<^j} to the desired accuracy. By in- 
tegrating Eq. (j3]) as shown in the Appendix, we get our 
main result 



(8) 



-cTx +^(t)cr^(a + a'^), (1) 



where the Pauli matrices are cix = \e){g\ + \g){e\^(jy = 
-i{\e){g\ - \g){e\),a^ = \e){e\ - \g){g\, while a is the an- 
nihilation operator for the harmonic oscillator. In a later 
section we will come back to some possible experimen- 
tal realizations of such time-varying coupling. For the 
purposes of this paper, we will consider the parametric 
coupling regime ^ = 0. In an interaction picture with 
respect to the free Hamiltonian, we have 



Hi = g{t)(Tz {ae 



-iQt 



(2) 



The open dynamics of the system is described by the 
master equation 

P = -i[Hi{t), p] + Cmp + Cqp + Cpp, (3) 

where the non-unitary contributions to the dynamics are 
given by the Lindblad operators 

^mP = ^(A^m + l)V[a]p^^NmV[o)]p, (4) 

= y(^. + ^m^-]P + ^NqV[cj^]p, (5) 



(6) 



In the above, cr_ = |^)(e|, cr+ = |e)(^|, V[A]p = 
2ApA^ — A^Ap — pA^A, where A is a generic operator, k 
is the damping rate of the oscillator, Fi and F2 are re- 
spectively the damping and dephasing rates of the qubit, 
while Nm (Nq) is the number of thermal excitations of 
the oscillator's (qubit's) environment. 



CHARACTERISTIC FUNCTION 
RECONSTRUCTION 

To begin our state reconstruction protocol, we assume 
to be able to initialize the system in the separable state 



where xiP) = {poD{P)} is the characteristic function. 



D{/3) 



being the displacement operator 



The notation {g^t) indicates that ^ and / are dependent 
on the specific realization of the time dependent coupling 
g{t)^ i.e. they are functionals. Their explicit form is 



a9.t) = 2i f dsg{s)e''''-^' (9) 

f{g,t) = jt^A{l-e--')\p{g,tf^KA [ ds\p{g,s)f 

Jo 

(10) 

/ dsg{s)e'^' sinh^s, (11) 

^0 



Ptot(O) = |+)(+| ^Po, 



(7) 



with 7 = Ti{Nq + 1/2) + 2F2 and A = Nm ^ 1/2. To 
reconstruct the characteristic function of the state po, 
one needs to "invert" Eq. (|9]), i.e. establish a mapping 
that associates any phase space point /3 to an appro- 
priate coupling and interaction time {g^^t^)^ such that 
i{gf3^tf3) = j3. One explicit example of such mapping is 
given in the next section. 

As shown in Eq. (jS]), the quantity that we can directly 
measure is given by the characteristic function evaluated 
at the phase space point ^(^, t) and "damped" by a factor 
^-fi9,t) (Jue to decoherence. Remarkably, this is still a 
valid representation of the state po, since the functional 
/(^, t) does not depend on the state of the oscillator and 
it is fully known in our theory once the coupling g(t) is 
assigned [21| . The actual value of the characteristic func- 
tion can thus be recovered even in the presence of finite 
decoherence. This however requires that we measure the 
expectation values in ([8]) with an accuracy greater than 
e"-^. Assuming that the relative error on (aj) scales as 
M~^/^, where M is the number of experimental runs, 
M ^ e^-^ is needed to measure the value of the charac- 
teristic function at the point ^ with sufficient accuracy. 
As we will see, / tends to increase as |^| is increased, 
meaning that decoherence imposes practical limits to our 
"reach" in phase space. 
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HARMONIC COUPLING 

An essential step of the reconstruction protocol is the 
inversion of Eq. (j9]). The realization of such inversion 
is clearly not unique and the specific form of g{t) and 
the choice of the interaction time can be optimized to 
best suit the available experimental apparatus. To give 
a concrete example, we will consider a coupling constant 
that oscillates harmonically at the mechanical frequency: 



2^' 



[ro + r sin ((/) - m)] , (12) 



where the exponential factor e'^^/^ is included to simplify 
calculations, r > 0, while ro is a constant that might be 
needed to keep ^r,(/)(^) inside the experimentally allowed 
range (for example, the constraint gr,cf){t) > would im- 
pose ro > r). For simplicity, we restrict the possible 
interaction times to integer multiples of the mechanical 
period. This choice is not mandatory, however it consid- 
erably simplifies the task of inverting Eq. (|9]), allowing 
us to keep the treatment analytical and compact. Eval- 
uating (j9]) for tn = n27r/l^, with n integer, gives 



(13) 



We see that the amplitude and phase of ^ are now related 
in a very simple way to the amplitude and phase of g(t). 
Had we chosen a generic interaction time, the value of 
would have depended in a more complicated way on the 
parameters {ro^r^ (j)^t). For high frequency nanomechan- 
ical oscillators it is often the case that ^max <C 1^ [ol, [iq| , 
thus the maximum value rmax of the parameter r will be 
typically small. A number of mechanical periods n > 1 
could then be used to reach phase space points of modu- 
lus greater than rmax- A possible inversion of Eq. ([9]) is 
then realized by associating the desired value of ^ to the 
corresponding parameters (r, 0, n), where n is the integer 
verifying (n - l)rmax < l^^l < ^^max, while r = \^\/n and 
= arg^. In practice the maximum number of cycles 
n and thus the modulus of ^ will be limited due to de- 
coherence, as we now show by explicit calculation. The 
damping exponent / defined in ([TO]) is in this case a com- 
plicated function of (n, r, ro, (/>). It is possible to obtain 
a compact result, valid for high quality factor oscillators, 
by approximating / to first order in hc/Q: 



f {gr,(ly, tn) ^ jtn + K^Atn 
^2 



sm ( 



+r^ 



ror- 



1 — 2n7r sin ( 



2n7r cos ( 



27r2 



cos ( 



47r2 



(14) 



If the parameters (ro, r, (j)) are kept fixed, it can be shown 
that (p!4j) is a strictly increasing function of the number 
of periods n, confirming the intuition that decoherence 
tends to be more severe for regions of phase space away 
from the origin [22[. As an example, if r^iax — 0.5, kA 




-2 



FIG. 1: Plot (a): Characteristic function of the Fock state 
|5). Plot (b): Phase-space representation of the same state as 
directly measured from the qubit, according to Eq.®. Plot 
(c) : Value of e^-^ as a function of the phase space point /3, giv- 
ing a lower bound to the number M of required experimental 
runs. Plots (b) and (c) are obtained with the damping expo- 
nent of Eq. (|14p . using the parameters ro = 0, r < 0.5, Q = 
27r X lOOMHz, n = 27t x 50kHz£7 = A = 27r x IMHz (these 
parameters are derived from 10] in the implementations sec- 
tion of this paper). We can notice how the damping factor 
e~-^ affects more significantly the characteristic function away 
from the origin, resulting in a larger number of required ex- 
perimental runs. Discontinuities in plots (b) and (c) are due 
to the discrete nature of the number of periods n in Eq. (O. 
As emphasized in the text, the characteristic function of plot 
(a) can be fully recovered from the measured data of plot (b), 
provided M ^ e^-^ runs are performed at each point /3. 



O.Oir^, 7 ~ O.Oir^, we can have n < 7 while still keeping 
the number of experimental runs reasonably low, since 
e'^f < 100 (recall that a number of experimental runs 
M e^^ is required to filter out decoherence from the 
measured data). In a later section of this paper we will 
show how this parameter range is easily achievable for 
the physical system proposed in [IQ]. This would allow to 
implement our protocol for phase space points of modulus 
1^1 = ^ 3.5, as we show in Fig. [T] for the particular 
case of the Fock state |5). 



MEASURING LOW ORDER MOMENTS 

It is often the case that full state reconstruction is not 
possible and only the characteristic function in the vicin- 
ity of the origin is availabe. This might still be sufficient 
to evaluate low order moments of the quadrature oper- 
ators Xe = ae~^^ ^ a^e^^ . To show this, we expand the 
right hand side of (|8|) in powers of r = |^|. Separating 
real and imaginary parts, we have a direct connection 
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between the oscillator moments and the expectations of 
pauli operators: 



(cT,(t))e/(3'*) = l 



1 



{ay{t))e 



f(g,t) _ 



riXe 



6 



1 

24' 



XI 



(15) 
(16) 



where 6 = arg{^} + 7r/2. Low order moments of an arbi- 
trary quadrature can therefore be obtained by first cor- 
recting the data for decoherence and subsequently per- 
forming a polynomial fit with respect to the variable r. 
The number of moments that can be reliably estimated 
with this method will depend on the available range of the 
parameter r. Useful information can be extracted even 
from the first few moments: while second order moments 
are sufficient to test squeezing, third and higher order 
moments can be used to investigate non-Gaussianity. 

SUPERPOSITION STATES 

One of the goals in experiments with nanomechani- 
cal oscillators is the preparation and verification of co- 
herent superpositions of classically distinct states, in or- 
der to explore the validity of the superposition princi- 
ple for macroscopic objects. In the previous sections, 
we fully addressed the problem of verification for arbi- 
trary states. Here, we show that the preparation of mo- 
tional superposition states can be achieved in the same 
framework. These possibilities together make our sys- 
tem a powerful toolbox for the investigation of nonclas- 
sical states of motion. We can prepare a superposition 
state in our system by following a few simple steps, sim- 
ilarly to what has been proposed in [g] and [13]. We 
suppose that the system is initialized in the state ([7j), 
and we assume that the oscillator has been pre-cooled 
to the ground state po = |0)(0|. At t = the cooling 
mechanism is switched off so that the time evolution of 
the system is described by ([3]). We let the coupled sys- 
tem evolve for a time t, then we measure the qubit in 
the basis \(f±) = ^ (|^) ± e^'^le)). In the absence of 
decoherence, the oscillator would be projected in the su- 
perposition state 



|V^Ja|a)±e-^'^|-a), 



(17) 



where | ±a) are coherent states with a = i dsg{s)e^^^ . 
The characteristic function is 



e 2 



cos(2^m{a/3*}) 



COS LP 



^- ^ • (18) 

2 + 2e-2|«Pcosv9 ^ ^ 

We can see that the diagonal terms | ± a){±a\ corre- 
spond to a sinusoidally modulated Gaussian peak cen- 
tered in zero, while the interference terms |±a)(=Fa| yield 
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FIG. 2: Comparison of the characteristic functions relative to 
the ideal superposition state and the prepared state. Plots (a) 
and (b) show respectively the real and imaginary parts of x for 
the superposition state |'0^,c/,) with a = 1, = |-. In plots (c) 
and (d) the same quantities are shown for the state prepared 
in the presence of decoherence, using the harmonic coupling 
of Eq.([T2|) and the parameters ro = 0, r = 0.5, n = 4, 11 = 
27T X lOOMHz, K = 27v X 50kHz, 7 = At A = 27r x IMHz. With 
the choice ip = 7r/2, the in-diagonal and off-diagonal terms of 
the density matrix are represented separately in the real and 
imaginary parts of x- Plot (c) shows that the state prepared 
with our method reproduces well the incoherent features of 
the superposition state, even though a slight steepening of 
the central peak due to heating can be seen in comparison to 
plot (a). On the other hand, the interference terms are more 
significantly damaged by decoherence, as it can be seen by 
comparing the imaginary parts (b) and (d). 



two Gaussian peaks centered at /3 = ±2a. The complex 
phase of these peaks is controlled by the relative phase cp. 
If we take into account finite decoherence, the described 
protocol yields an imperfect superposition state (related 
calculational details can be found in the Appendix). In 
particular, we can expect a steepening of the central peak 
due to heating (i.e. the state becomes mixed) as well as 
a reduction in the height of the interference peaks, due 
to dephasing. This is shown in Fig. [21 where the charac- 
teristic functions obtained with our protocol in the ideal 
and finite-decoherence cases are compared. To conclude 
the present section, we point out the dual behavior of 
the characteristic function as compared to the Wigner 
function [sf. In the latter representation, interference 
terms appear as oscillations, while diagonal terms yield 
non-centered Gaussian peaks. Moreover, in the Wigner 
function heating results in broadening of the Gaussian 
peaks rather than steepening. These observations are all 
consistent with the fact that the two representations are 
connected by a symplectic Fourier transform. 
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POSSIBLE IMPLEMENTATIONS 

Being rather general, our reconstruction method can 
be apphed to any experimental setup in which the qubit- 
oscillator coupling can be coherently tuned over time. 
We shall give some examples of realistic settings in which 
such level of control can be achieved. 

As our main example, we consider the theoretical pro- 
)osal by Fei Xue et at. in the field of nanomechanics 
IQl . In their paper, the coupling between a flux qubit 
16| and a nanomechanical resonator is controlled via the 
amplitude of an external magnetic field, according to 



g{t) = rjB{t), 



(19) 



where is a constant that depends on the specific sys- 
tem realization, while B(t) is the magnetic field am- 
plitude along an appropriate direction. They estimate 
?7 27r • 0.8MHz/mT for a realistic choice of parame- 
ters, emphasizing how the control of such coupling can 
be achieved without interfering with the non-interacting 
part of the system Hamiltonian. From the form of Eq. 
(p!9|) . we see that both the magnitude and sign of the 
coupling can be controlled over time by tuning the ex- 
ternal magnetic field. Still referring to [lo[, we assume 
that the oscillator has a frequency 17 :^ 27r x lOOMHz, 
ground state spread xq — 2.6 x 10~^^m and quality fac- 
tor Q — 2 X 10^, resulting in 27r x 50kHz and 
HiA 27r X IMHz at a temperature T ~ lOOmK. Assum- 
ing a maximum magnetic field intensity -B^ax 

lOmT, 

we get ^max — 27r X 8MHz. With such parameters, 
in ([12]) we can set ro = and r < 27r^max/^ — 0.5. 
The implementation of the harmonic time dependence of 
([12]) requires manipulation of the magnetic field source 
currents at radio-frequency, which can be achieved with 
modest technology. For the flux qubit we take Fi = F2 = 
27r X 0.4MHz and TVg :^ 0, yielding a total dephasing rate 
7 = Fi/2 + 2F2 27r X IMHz. Inserting these parame- 
ters in Eqs. (fT3|) and ([T4|) . considering rmax — 0.5, n = 7, 
and assuming the Markovian approximation to be valid, 
we are able to reach phase space points up to a distance 
1^1 3.5 with a relatively low number of experimental 
runs, since e^-^ < 100 as shown in Fig. [TJ This parame- 
ter range is already sufficient to identify a large variety 
of states, including Fock states of low order, coherent 
states and superposition states with |a| < 1. To measure 
states with a thermal excitation, even smaller values of 
1^1 could be sufficient, since finite temperature induces 
localization of the characteristic function around the ori- 
gin. As the number of mechanical periods is increased, 
further regions of phase space become available, however 
the number of required runs blows up rapidly for n > 7, 
and already for n = 10 we have e^-^ - 10^ With the 
above parameters it is also possible to prepare superpo- 
sition states with \a\ 1, as shown in Fig. O 

As a second example we consider the recent experi- 
mental work of Srinivasan et al. [l2|. In it, the au- 



thors demonstrate a novel circuit-QED architecture, in 
which both the internal levels splitting of a qubit and 
its dipolar coupling to a microwave cavity mode can be 
independently tuned over a wide range of parameters. 
The system was theoretically proposed in [11|. By con- 
trolling external bias voltages, the authors are able to 
continuously vary g in the range 2it • 200kHz^46MHz. 
The harmonic oscillator is in this case a microwave cav- 
ity mode of frequency 2it - 5GHz, so that the ratio 
g/Q. varies between a negligible value and ~ 0.01. Re- 
cent results from the same group suggest the possibility 
to reach g ^ 27t ■ 300MHz, which would push the ratio 
up to ~ 0.06. Even though the experiment by Srini- 
vasan focuses on the Jaynes-Cummings regime, we argue 
that there is no fundamental reason why the same ideas 
should not work in the dispersive regime required by our 
protocol, and we hope to see experimental confirmations 
in this direction in the near future. 

The application of our scheme to systems based on 
charge qubits ^ is more problematic at the actual state 
of technology, essentially due to fast qubit dephasing. 
Moreover, the background charge noise acting on the 
qubit is usually not well described by a Markovian master 
equation of the form ([3]). However, the situation might 
be improved by combining our protocol with charge-echo 
techniques jl7|, which we leave for future investigations. 



CONCLUSIONS 

Before concluding, we emphasize that estimating deco- 
herence through Eq. ([3]) is in general an approximation, 
since many environments are not exactly Markovian. For 
environments that show non-negligible deviations from 
Markovianity, we can expect the predictions of our model 
to get progressively worse as the interaction time is in- 
creased. To correct this, non Markovian effects could be 
included numerically in the model. 

Even when decoherence is accurately modeled by the 
master equation ([3]), in a real experiment unavoidable er- 
rors will arise due to random fluctuations in the coupling 
strength g and a limited accuracy in the control of the 
interaction time t. Similarly, the decoherence parameters 
appearing in the master equation might be known with 
a non- negligible uncertainty. However, such errors can 
be bounded, due to the continuity of the functionals in 
Eqs. ([9l fTT]) with respect to (g^t) and the decoherence 
parameters. 

To summarize, we presented a scheme in which a para- 
metrically coupled qubit can be used to measure the char- 
acteristic function of a nanomechanical oscillator. By 
introducing the possibility of a time-varying coupling, 
we have shown how the characteristc function can be 
measured just by exploiting the system's time evolution, 
while the effect of Markovian decoherence can be filtered 
out by increasing the number of experimental runs. 
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APPENDIX: SOLUTION METHODS 

We solve the master equation Q b y using phase space 
methods similar to those used in [18| . We consider a rep- 
resentation in which a matricial characteristic function 
is used to describe the state of the coupled system. We 
decompose the total density matrix at time t as 

Ptotit) = Pe{t) \e){e\ + pg{t) \g){g\^ 

-^P^{t)0\e){g\^p_{t)0\g){e\. (20) 

By defining the characteristic function for each element 
as Xj(/3, = tra {pj{t)D{f3)}^ we can define the matricial 
characteristic function as: 



Xtot {P,t)= Xe (/?, t)\e){e\+ Xg {P,t)\g){g\ + 
+ X+W,t)\e){g\+x-W,t)\g){e\. 



(21) 



At this point, we have to convert ([3]) to a system of cou- 
pled partial differential equations for the functions Xj^ 
which can be done by using standard techniques [sl . 

The expectation values required for our state recon- 
struction protocol are: 

{a,{t)) = tr{ptot(tK} = X+(0,t) + x-(0,t), (22) 
{ay{t)) = tT{ptot{t)cTy} = i [x+(0,t) - x-(0,t)] , (23) 

therefore for the time being we only need to compute 
the evolution of the off diagonal elements x±(/^5^)- The 
equations for x± are already in diagonal form: 



dtx± = ±2i5(t)(e-'"*a^. - e^"*a^)x± + C^X± - 7X±, 

(24) 

where the differential form of the mechanical Lindblad 
operator is CmXj = "f {f^dp + P*dp, + 2A|/?|2) Xj. The 
solution of Eq. is 



X±(/3, t) = x±(/3e- 1* T o)e-^(i-^""')l''T''l'-^ (25) 

where v = ^ k,A ds\p{g^ s)^, and ^, p are the func- 
tionals defined in (0) and ([II]). Now, Eq. (jH]) follows by 
considering the initial conditions corresponding to the 
initial state (O, i.e. x+(/5,0) = X-(/5,0) = |x(/5)- Note 
that the factorization of the effect of decoherence in the 
right hand side of Eq. (|8]) is in strong analogy with 
the beam-splitter model for the decoherence of a bosonic 
mode [l9| . 



To treat the superposition state preparation, we begin 
by calculating the characteristic function at time t after 
postselection of the outcome \^±) 

{V±\XtotW,t)\ip±) = XeW,t) + Xg{P,t) + 

±e-'^X+(/?,t) + ±e'^X-(Ai). (26) 

After normalization, we get the characteristic function of 
the imperfect superposition state 



X± 



Xe{P,t)+XgiP,t)±e-^'^X+iP,t) ± e"^X-iP,t) 
2±e-^Vx+{0,t)±e'^X-{0,t) 



(27) 

We see from the above equations that we need to com- 
pute the evolution of the remaining elements Xe,g- For 
simplicity, we neglect thermal excitations in the qubit 
bath and take Nq 0. This is justified by the fact that 
the qubit frequency can be several orders higher than 
the oscillator's. With this approximation, the equations 
of motion are 



dtXg 



ftXe = ^^(t)(e-^^^/3 + e^^^/3*)xe 



(28) 

CmXe-^lXe- (29) 

Eq. ([29|) is homogeneous and can be readily solved 

Xe(/3,t) = e-^^*-^0--"*)l^l^+^'''-^*^e(/3e-t*,0), 

(30) 

where A = ie~^'^ dsg{s)e^^^~^^^ (note that this is also a 
functional depending on g and t). Eq. ([28|) is non homo- 
geneous due to the presence of the term FiXe, therefore 
yielding a more involved solution 



Xg{P,t)=Xg{P,t) 



l + Ti 



*^^^ Xe(/3et(^-'),s) 
''xg(/?et(--*),s) 



, (31) 



where Xg is the homogeneous solution of Eq. (f28|) . 

XgiM = e-^i'-^-^')\^\'+''^-'^\,{Pe-^\0). (32) 

Since we supposed that the oscillator is initialized in 
the ground state, the initial conditions are now given 
by x,(/3,0) = |e-*l^l'forj = e,^,+,-. 
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